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nomenon of power-law level repulsion ||19| , |2T|j which is globally very well described 
by the Brody distribution. 
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The study of energy level statistics of generic quantum Hamiltonian systems 
whose classical dynamics is between integrability and full chaos persists to be a 



challenging problem [19, 20, 21]. Until very recently there were incompatible results 
concerning the so called level spacing distribution P{S) where P(S)dS is a probabil- 
ity that a randomly chosen spacing between two adjacent energy levels lies between 
S — dS/2 and S + dS/2. Berry and Robnik || derived the semiclassical level spacing 
distribution P BR (S) assuming the principle of uniform semiclassical condensation 



[g, [23|, 16, 11, of eigenstates onto classical invariant components (which can 
be either regular - tori or irregular - chaotic) and the statistical independence of 
the level subsequences belonging to various disjoint classical invariant components. 
(Regular levels associated with quantized invariant tori may be merged together 
giving the well known Poisson distribution Pp i S son('S') = e~ s .) Using the factoriza- 
tion of gap distributions E(S) = do{p — S)P(S) upon statistically independent 
superposition of spectra one may write 

E™(S) = E Poisson (p 1 S)E GOE (p 2 S) (1) 

for the simplest case of only one chaotic component with relative measure p 2 and 
regular components with total relative measure p\ = \ — p%. The Berry- Robnik 
distribution does not exhibit level repulsion, since -P^ R (0) = 1— p\ ^ 0. On the other 
hand there has been a vast amount of phenomenological evidence [^] in favour of 
the so called fractional power law level repulsion which is globally very well described 
by the Brody |J distribution 

p£(S) = aS?eM-bS l3+1 ), a = (P + l)b, b = [T(l + (0 + l)" 1 )]^ 1 (2) 

or the more sophisticated Izrailev |TTJ distribution which are characterized by the 
noninteger exponent (3, P(S — > 0) oc S?. Numerical spectra which contain even 
up to several ten thousands energy levels of quantum Hamiltonian systems with 
mixed classical dynamics typically still exhibit the phenomenon of fractional level 
repulsion, with statistically significant global fit by the Brody distribution. In such 
cases there was a persisting puzzle how the level spacing distribution converges to 
the semiclassical Berry-Robnik distribution as one increases the sequential quantum 
number or decreases the value of effective H. However, recently we have succeeded 
to demonstrate the ultimate semiclassical Berry-Robnik level spacing distribution 
in a rather abstract 1-dim time-dependent dynamical system, namely the standard 
map on a torus, and showed the smooth transition from Brody-like to Berry-Robnik 
distribution as Ti decreases |H| ^] (see also |13|0 . The transition was excellently 
described by two parameter (pi,(3) Berry- Robnik-Brody model in which we substi- 
tute the GOE model for the chaotic part E GOE (p2S) in Berry-Robnik formula ([!]) 
by the Brody model E R (p 2 S) with some exponent (3. The major goal of this letter 
is to demonstrate this scenario in a generic 2-dim autonomous Hamiltonian system. 

Another goal of this letter is to demonstrate the practical power of the recently 
developed quantum surface of section method [fTS], [TJ], ^| (which has been motivated 
by the semiclassical version developed in [ffl and extensively numerically investigated 
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in H) and whose most thorough and complete presentation so far is given in |15| . We 
shall use the reactance matrix formulation of the quantization condition which has 
practical advantages over the scattering matrix formulation in the case of semisep- 
arable systems ( ||15| ) section 2.7). Here we give a brief and hence rather heuristic 
description of the method for the quantum Hamiltonians H in 2-dimensional con- 
figuration space (CS) with coordinates (x, y) where the line y = represents the 
configurational surface of section (CSOS) while we use Dirichlet boundary condi- 
tions on the boundary lines y = > and y = —b± < 0. Let ^ an (x,y,E) be the 
solutions of the Schrodinger equation H^ an (x,y, E) = E^ an (x,y,E) on the upper 
(y > 0, a =t= +) / lower (y < 0, o =|= — ) side of CS which satisfy the boundary 
conditions ^/ an (x, 0, E) = u n (x), ^/ an (x, ab a , E) = 0, and u n (x) is some complete 
set of functions for the small Hilbert space of L 2 functions over 1-dim CSOS, e.g 
the eigenfunction of the reduced Hamiltonian H' = H\ y=0 , H'u n (x) = E' n u n (x). 
The full eigenfunction ty(x,y,E) of H can be expanded in terms of partial eigen- 
functions ^ an (x, y, E) on both sides as ^(x,y,E) = I] n c - n |fc n (-E')|~ 1//2 \l/ - n (x, E), 
where a = sign(y) while the square roots of the wavenumbers k n = h^ 1 ^ 2m(E — E' n ) 
provide a useful normalization. In order that ^(x, y, E) would be a nontrivial eigen- 
function on the entire CS it should be continuously differentiable on CSOS (at 
y = 0). Using the completeness of the set u n (x) the requirement for continuity 
gives Cf n = C| n , whereas requiring continuity of the normal derivative yields the 
singularity condition for the real symmetric reactance matrix R = Rj + Rj_, 

R an i{E) = a\k n {E)k l {E)\- 1 ' 2 j ' dxV an {x,0)d y * al {x,0), (3) 
detR(E) = 0. (4) 

This equation is equivalent to the more physical but numerically less effective (due to 
complex non-symmetric arithmetic) quantization condition det(l — T±(E)T^(E)) = 
|T2|, [TJ], [23], |T3| in terms of generalized (non-unitary) scattering matrices T a (E) 



of the two scattering problems (obtained by cutting one half of the CS along CSOS 
off and attaching the waveguide (flat in the y-direction) instead) which have a finite 
number N a of propagating - open, and infinitely many evanescent - closed modes 
e ±tkn ^ y u n (x), for k 2 (E) > 0, and k^E) < 0, respectively. The scattering matrices 
are related to [non-real) reactance matrices, by = (1 + zR cr )(l — iR,,) -1 , where 
the latter are made real by a simple diagonal transformation Ro- = DR CT D, where 
D = diag(l, 1 . . . N times ... 1, y/i, y/i . . .). 

We have applied this method to a semiseparable system which is separable 
above/below CSOS but not separable on the whole CS, namely to 2-dim semisepa- 
rable oscillator (SSO) with the Hamiltonian 

H = -\h 2 (d 2 y + dl) + \(x-\s\gn(y)a)\ -6 X < y < b h (5) 

with the parameters a, b a , Ti. SSO has a scaling symmetry (a, b a , h, E) — > (ao, ab a , a 2 h, a 2 E). 
The reduced Hamiltonian is just a simple 1-dim harmonic oscillator —\h 2 d 2 + \x 2 
with eigenf unctions u n (x) = (\pKfil n w)?)~ l l 2 exp(— x 2 /2h)H n {x/y/h) and eigenener- 



gies E' n = (n+7;)h determining the wavenumbers k n (E) = h 1 \/2E — (2n + l)h, n 
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0, 1 ... It is easy to derive an explicit expression for the reactance matrices for SSO 
R T (£) = J(£)OC t (£)O t J(£), R 4 (£) = J(E)0 T C i ( J E;)OJ( J E;), (6) 
where J(E),C a (E) are rea/ diagonal matrices 

3 nl (E) = 5 nl \k n (E)\-^ 2 , C anl {E) = 8 nl k n (E) cot(k n (E)b a ) , (7) 
and O is real orthogonal shift matrix 

O n i = J dxu n (x)ui(x + \a) (8) 
whose matrix elements can be calculated via numerically stable symmetric recursion 

- . - . „ . -IV 




2^/ \nj ' VS2hn V32M 

It is important how to truncate these infinitely dimensional matrices for the numer- 
ical calculation. One has to consider all the N = io\md(E/h) open modes plus as 
many N c closed modes so that the numerical results (roots of eq. (||)) converge. I 
have used semiclassical arguments (the SOS (x,p x ) phase space supports of coherent 
state representation of the states u n (x), n = 1, . . . , N Q + N C should cover the supports 
of the states ui(x + |a), I = 1, . . . , N Q ) to estimate the minimal number of closed 
modes 

/2a a 2 \ 
N c w -= + — JV . (9) 



2E 2E t 

Dimension of matrices = N a + N c is thus usually (for small a) only little larger 
than the number of open modes N Q . 

It is also very important to stress that the shift matrix and therefore also the 
reactance matrices are effectively banded. I have obtained semiclassical formula for 
their bandwidths using overlap condition for the coherent state representation of the 
SOS-states u n (x) and U\{x + |a) 

bandwidth(R (T ) = 2 bandwidth(O) w %\FlE. (10) 

n 

Note that the function f(E) = detR(^) has singularities (poles) at the points E 
where for some n, k n (E)b a is a multiple of ir. But between the two successive poles 
f(E) is smooth (even analytic) real function of real energy E. I have devised an 
algorithm for calculation of almost all levels — zeros of f(E) within a given interval 
[Ei,Ef] which needs to evaluate f(E) only about 25 times per mean level spacing 
while it typically misses less than 0.5% of all levels. The control over missed levels 
is in general very difficult problem. The number of all energy levels below a given 
energy E, M{E) can be estimated by means of the Thomas-Fermi rule 

ME) w M F (E) = h -±±^(2Ef 2 = 0(N 2 ). (11) 

3nh 
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But this formula is generally not very helpful even if next semiclassical corrections 
are negligible since the fluctuation of the number of levels in an interval [Ei, Ef\ 
is proportional to Jl\f{Ef) — M{Ei) except in the extreme case of fully chaotic 
systems where the spectra are much stiffer and the fluctuation is proportional to 
\og[H{Ef) — Af(Ei)] so that Thomas-Fermi rule can be used to detect even single 
missing level 0, |T|. 



We have chosen the following values of the parameters for our numerical demon- 
stration a = 0.03, &| = 5.0, b[ = 10.0, E = 0.5 while for quantal calculations we 
take energy to be in narrow interval around E = 0.5. For illustration we plot the 
classical SOS (x,p x ,y = 0) in figure 1. There is only one dominating chaotic com- 
ponent with relative measure p2 = 0.709 ± 0.001 and the regular region still with 
some very small chaotic components with the complementary total relative measure 
Pi = 1 — p 2 — 0.291 ± 0.001. For the quantal calculations we have chosen two 
different values of h = 0.01 and h = 0.0003 which correspond to sequential numbers 
M « 16 000 and Af « 17684 000, respectively. 

In the first case (h = 0.01, M ~ 1.6 ■ 10 4 ) we have calculated 14 231 levels in 
the interval 0.35 < E < 0.65. We have performed a y 2 test and obtained statisti- 
cally significant fit of P(S) by the Brody distribution with (3 = 0.142, Xb = 5320 
and nonsignificant Berry-Robnik fit with p\ = 0.548, Xbr = 130 000 (see figure 
2). In order to present most detailed information we plot cumulative level spac- 



ing distribution W(S) = f$ daP(a) and deviation of numerical U— function [IS 



U(W(S)) = (2/7r)arccos(y l-W(S)) from the best fit Berry-Robnik U-function 
U(W BR (S)) versus W(S) which has a nice property that the estimated statistical 
error SU = 1/tc\/N and the density of numerical points along abscissa are constant. 
In spite of the already very high sequential number this is still an example of the 
so-called near-semiclassical regime characterized by the fractional power law level 
repulsion. 

In the second case (h = 0.0003, M ~ 1.8 ■ 10 7 ) we have calculated 13 445 levels 
on the interval 0.49985 < E < 0.500105 and found significant fit by the semiclas- 
sical Berry-Robnik formula (see figure 3) with the correct value of regular volume 
pi = 0.286 ± 0.005, Xbr = 12 150 while Brody fit becomes highly nonsignificant, 
(3 = 0.367, = 249 000. Thus we have demonstrated the so-called far semiclas- 
sical regime with the quantum value of p\ which excellently agrees with the clas- 
sical regular volume (the small deviation is within error bars). Fit to combined 
Berry- Robnik-Brody model does not not significantly improve x 2 = 11950 while it 
substitutes GOE model for the chaotic part by the Brody model with f3 ~ 0.85. 

Large square root number fluctuations prevent to determine the number of missed 
levels by using Thomas-Fermi rule (although higher order semiclassical corrections 
are negligible in this regime). One can compare the number of levels M{E) with the 
number of levels Nq{E) or Afoo(E) for the two nearby integrable - separable cases 
(with the same b a but with a = (single box limit) or a — > oo (two box limit), 
respectively) since the leading order semiclassics (Thomas-Fermi rule) does not de- 
pend upon defect a. Mq{E) and Noo(E) can be easily calculated numerically and 
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large scale fluctuations of Af(E) — A/'o i oo(-E') turn out to be much smaller than the 
fluctuations of J\f{E) — Af TF (E) suggesting that we have missed less than 20 levels 
out of 14 231 at % = 0.01 (figure 2) and 40 - 80 levels out of 13 445 at h = 0.0003 
(figure 3). Note that in the first case {Ti = 0.01) there was much less almost degen- 
erate pairs of levels (and therefore missed levels) due to the level repulsion. 



In conclusion I should emphasize that the present work presents clear demonstration 
of the semiclassical Berry-Robnik level spacing distribution in a generic 2-dim au- 
tonomous Hamiltonian system between integrability and chaos, i.e. semiseparable 
oscillator. This would not be possible without application of the quantum surface 
of section method which enabled us to calculate 13 445 consecutive levels with se- 
quential numbers M « 1.8 • 10 7 within a week of Convex 3680 CPU time. For 
this system, quantum surface of section method requires 0(a 2 J\f 3 ^ 2 ) FPO/level. In 
the forthcoming publications I will |DJ discuss in detail the energy spectra, quan- 
tum eigenstates, and quantum SOS evolution by the quantum SOS method in the 
semiseparable oscillator, and [T?| apply the quantum SOS method to a more realistic 
example, namely the diamagnetic Kepler problem |fj, [7], |22]| . 
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Figures 



Figure 1 Classical surface of section y = with coordinates x, p x of the semisepara- 
ble oscillator a = 0.03, 6 T = 5.0, b l = 10.0, E = 0.5. 250 orbits with 40 000 crossings 
of SOS each are shown. 

Figure 2 Cumulative level spacing distribution W(S) (a) and deviation of its U- 
function from the best fit Berry-Robnik distribution (b) for 14 231 consecutive levels 
of SSO at a = 0.03, 6 T = 5.0, b i = 10.0, 0.35 < E < 0.65, h = 0.01. Thick full curves 
(a,b) (within ± one sigma error band (b)) represent numerical data, the thin full 
curve (a) is best fit Berry-Robnik distribution with pi = 0.548 whereas the dashed 
curves (a,b) is the best fit Brody distribution with (3 = 0.142. The dotted curves 
(a) represent the limiting Poisson and GOE distributions whereas the dash-dotted 
curves (b) represent the nearby Berry-Robnik curves with p 1 = 0.548 ±0.01, ±0.02. 

Figure 3 The same as in figure 2 but now at h — 0.0003 for 13 445 levels in the 
interval 0.49985 < E < 0.500105 living in the true (far) semiclassical regime (see 
text). The Berry-Robnik distribution with pi = 0.287 is now statistically significant 
and for illustration of accuracy of the fitted p\ we provide also Berry-Robnik curves 
for pi = 0.287 ± 0.01, ±0.02 (dash-dotted curves (b)). 
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